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We propose variational states for the ground state and the low-energy collective rotator excitations 
in negatively charged C^>~ ions (N = 1 ... 5). The approach includes the linear electron-phonon 
coupling and the Coulomb interaction on the same level. The electron-phonon coupling is treated 
within the effective mode approximation (EMA) which yields the linear ti u ®H g Jahn- Teller problem 
whereas the Coulomb interaction gives rise to Hund's rule coupling for N — 2, 3, 4. The Hamiltonian 
has accidental SO (3) symmetry which allows an elegant formulation in terms of angular momenta. 
Trial states are constructed from coherent states and using projection operators onto angular mo- 
mentum subspaces which results in good variational states for the complete parameter range. The 
evaluation of the corresponding energies is to a large extent analytical. We use the approach for 
a detailed analysis of the competition between Jahn- Teller effect and Hund's rule coupling, which 
determines the spin state for N = 2,3,4. We calculate the low-spin/high-spin gap for N = 2, 3,4 
as a function of the Hund's rule coupling constant J. We find that the experimentally measured 
gaps suggest a coupling constant in the range J = 60 — 80 meV. Using a finite value for J, we 
recalculate the ground state energies of the C^ - ions and find that the Jahn- Teller energy gain is 
partly counterbalanced by the Hund's rule coupling. In particular, the ground state energies for 
TV = 2, 3, 4 are almost equal. 



PACS numbers: 71.7D.Ej, 73.61.Wp 
I. INTRODUCTION 

The t\ u ® H g Jahn- Teller problem, where electrons 
in a threefold degenerate orbital interact with a five- 
fold degenerate phonon multiplet, is known since more 
than 30 years. It first arised for the particular case of 
p-electrons in a cubic systems which are equally cou- 
pled to Eg and T^g vibrational modesi. On the level 
of linear coupling, an equivalent problem arises in neg- 
atively charged C^, - ions (N = 1...5). These mate- 
rials experienced particular interest when superconduc- 
tivity was observed in alkali-doped A 3 Ceo (A=K,Cs,Rb, 
for a review see Ref . |U) ■ The neutral Ceo molecule is a 
closed shell system and highly symmetric. It has icosa- 
hedral symmetry which is the largest threedimensional 
point-group with 1-, 3-, 4- and 5-dimensional irreducible 
representations (IR) 3 . The lowest unoccupied molecu- 
lar orbital (LUMO) of Cgo is threefold degenerate and 
has t\ u symmetry. It couples to two non-degenerate A g 
phonon modes and eight 5-fold degenerate H g phonon 
multiplets^. In the present work we focus on the non- 
trivial coupling to the H g multiplets. We restrict our 
attention to linear coupling and approximate the eight 
H g multiplets by one effective multiplet which gives rise 
to the linear, single-mode t\ u (£> H g Jahn- Teller problem. 
Furthermore, we will use the fact that the icosahedral 
IR's tiu and H g correspond to the L — 1, 2 IR of SO (3) 
which don't split under the icosahedral symmetry 3 -. As 
a consequence, the linear t± u <g> H g Jahn- Teller problem 
is equivalent to the problem of p-electrons interacting 
with d-phonons and recovers accidental SO (3) symme- 
try. Therefore, the present work really treats the linear 
p®d Jahn- Teller problem and we will mostly speak of p- 
electrons and d-phonons throughout the following. Note 



that the accidental SO (3) symmetry would be lifted by 
higher order coupling terms^i In addition to electron- 
phonon coupling, electrons in Cgo also interact via the 
Coulomb interaction. This leads to the so-called Hund's 
rule coupling (see Ref. H and references therein). Below, 
we will consider both, the Jahn- Teller and Hund's rule 
interaction, and discuss the competition between them. 

The p <g> d Jahn- Teller cannot be trivially solved as, 
for example, a single displaced harmonic oscillator. This 
led to various different approaches. Early work on the 
case N = 1 was done by O'Brie n 1 ' 9 ' 10 . The first treat- 
ment for all fillings was carried out in two subsequent 
papers by Auerbach, Manini and Tosatt i 11 ' 12 . The first 
paper is based on a semiclassical (also called adiabatic) 
approximation which yields the effective Hamiltonian in 
the strong coupling limit^i. The intermediate regime is 
explored using exact diagonalization. The weak coupling 
limit is treated in the second paper using perturbation 
theory 12 . The effect of Hund's-rule coupling was studied 
subsequently using the same approach— 3 -. These works 
led to a good understanding of the t\ u <g> H g Jahn- Teller 
problem. A complete discussion is given in Ref. 0. 

However, none of the schemes just discussed applies 
to the whole coupling range. Even "exact" diagonaliza- 
tion is only valid for the small and intermediate regime 
because it suffers from truncation of the phonon Hilbert 
space at high coupling. Moreover, the approaches don't 
provide wavefunctions in the intermediate regime, which 
is precisely the regime of Cgo- There have been var- 
ious attempts to construct wavefunctions for all cou- 
pling regimes on the basis of coherent states and pro- 
jection technique o 14 ' 15 ' 16 ' 17 ' 18 . The formalisms used were 
rather involved and led to complicated analytical expres- 
sions and multidimensional integrals. Here, we propose 
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a mathematically equivalent but much more convenient 
formalism which is based on the use of projection opera- 
tors. 



We construct variational states in two steps: First, we 
start with a product state \^(q)) = \ip) e <8> \q) p which 
minimizes the electron-phonon coupling. The phonon 
part \q) p is a coherent state which corresponds to the 
displacement q of one phonon coordinate. This displace- 
ment leads to a splitting of the degenerate electronic lev- 
els. The electron part \ip) e is chosen such as to min- 
imize the energy of the electrons for a given splitting. 
The state \^>) is not SO (3) symmetric and also not a 
good angular momentum state. However, the SO(3) sym- 
metry of the Hamiltonian requires eigenstates to be an- 
gular momentum states. Therefore, we use projection 



II. HAMILTONIAN 



The Hamiltonian describing the full multi-mode Jahn- 
Teller problem with Coulomb interaction of the ion 
has four terms: 



H — H, p + H ep + Hj + H 



(1) 



The hrst term is the energy of the 8 d-phonon multiplets 



(2) 



where uj a are the frequencies of the phonon multiplet: 
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(a = 1 



and aka are the phonon construction 



and annihilation operators. According to the definition 
in appendix [XJ the quantum numbers k — — 2 . . . 2 corre- 
operators Q L MKl as defined in (23]), to construct angu- spond to the d-symmetries V3xy, y/3yz, z 2 - (x 2 + y 2 )/2, 

\fZxz, \^3(x 2 — y 2 )/2 respectively. As discussed in the in- 
troduction, we only consider linear electron-phonon cou- 
pling which is given by 



lar momentum states \LMK,q) = Qmk\^{i))- State 
\LMK, q) is the final variational wavefunction with good 
quantum numbers LM and one variational parameter q. 
The expectation value of H with respect to this state is 
{MQkmHQmk^) = ^\HQ L KK \il>). The equality arises 
because H is a scalar and commutes with the projection 
operator. This is an essential simplification because each 
projection operator carries an integration over Euler an- 
gles, as can be seen from definition (J23J. Previous works 
haven't made use of this property which, as we show be- 
low, allows an almost analytical treatment of the prob- 
lem. 



The present approach is interesting because it is to 
a large extent analytical and applies to the whole cou- 
pling range. Energies can be calculated with moderate 
effort which allows a detailed analysis of the competi- 
tion between Jahn- Teller effect and Hund's rule coupling. 
This competition determines the spin configuration of the 
ground state for the cases N = 2, 3, 4 which we calculate 
for the complete parameter range. 



The paper is organized as follows: In section [IT] we in- 
troduce the Hamiltonian as well as the effective mode 
approximation which reduces the multi-mode problem to 
a single-mode problem. The properties of the phonon co- 
herent states are discussed in section lHll Scction[TV]to[VT] 
treat the cases N = 1, 2, 3 respectively. As discussed be- 
low, the cases N — 4, 5 are equivalent to N — 2, 1 due to 
particle-hole symmetry. Results specific to parameters of 
Ceo are discussed in section IV*lT1 Due to the SO (3) sym- 
metry, angular momenta and its eigenstates, the spher- 
ical harmonics, play an important role in the present 
work. We always use real spherical harmonics Ylm(&, 4 1 ) 
which have cos(Mcf)) or sin(M<f>) dependence. Rotation 
and products of the real spherical harmonics lead then to 
real Wigner-D functions D\ lK and new Clebsch-Gordan 
coefficients denoted by fi^S M . The corresponding 
definitions are discussed in detail in Appendix lAl 
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qkc 



(3) 



where g a are the coupling constants and <jka — i a \ a + 
afea)/v / 2 the operators for the phonon coordinates. cJj S 
and c n ' s are the electron operators. They have spin s and 
quantum numbers n = —1,0, 1 which correspond to the 
p-symmetries y, z, x respectively. R 2t ^ ln , are the Clebsch- 
Gordan coefficients for the real spherical harmonics (see 
appendix [X| . The resulting matrix elements of H ep are 
given in (TTB]) . Note that both, the electron operators cJj S 
and c n ' s as well as the phonon operators a\. a , a,k a and 
qka are tensor operators with rank 1 and 2 respectively. 
The sum in involving the Clebsch-Gordan coefficients 
is the simplest non-trivial scalar which can be built from 
rank 1 and rank 2 tensors. 

The last two terms in (fTJ) describe the Coulomb inter- 
action among the p-electrons. The Coulomb interaction 
splits the charge states N = 2,3,4 into multiplets char- 
acterized by the total electron angular momentum and 
the spin because spin-orbit coupling is omitted. All mul- 
tiplet energies can be expressed by two parameters J and 
U giving rise to the two terms Hj and Hjj^ . They have 
the form 



J 



U 



H v = -jV(iV-l) 



y ] ( C ns C ns' Cms ' Cms C ns C ms' Cns ' Cms ) '(^) 



(5) 



where N is the total number of electrons. In the litera- 
ture, Hj is referred to as the Hund's rule coupling and 
leads to the multiplet splittings listed in Tab. HI The 
second parameter, U , is the overall charging energy. It 
is of the order [7=1 — 3 eV, depending on the screen- 
ing£i2£. J is at least one order of magnitude smaller 8 as 
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-J 


3 


3/2 
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TABLE I: Multiplet energies £ mu it of C^ ~ for JV = 2, 3, 4 
electrons which arise from the Hund's rule coupling Hj given 
in (j4| . The quantum numbers N, S, L £ denote the number of 
electrons, the total spin and total electron angular momen- 
tum. 



will be discussed in section IVTfl In the following we will 
drop Hjj because we always work with a fixed number of 
electrons. 

All terms in the Hamiltonian {1} are scalars and there- 
fore SO (3) invariant. Hence, the total angular momen- 
tum L = L e + L pQ , which is the sum of the electron 
and phonon angular momenta, is conserved and eigen- 
states of H have quantum numbers L and M. In ad- 
dition, H is particle-hole symmetric if Hjj is neglected. 
Therefore it is enough to study the cases N — 1, 2, 3. 

The Hamiltonian above can be simplified by introduc- 
ing one effective d-phonon multiplet instead of 8 multi- 
plets. This results in the effective mode approximation 
(EMA) which was in detail investigated by O'Brien 10 . 
The phonon operators for the effective mode are a super- 
position of the original modes: 

8 8 

ak = u " a ka with y]u% = l. (6) 

a— 1 a— 1 

The real coefficients u a can be determined by the vari- 
ational principle: If 1$) is a state which only contains 
excitations of the effective mode, then the following gen- 
eral relation holds (for a derivation see Ref . [2lh : 

mm(H) 9 =u ~u>+u(Hf+H^)s+(Hj+H u )s, (7) 

{u a } z 

where loq = 5/2 ^ Q uj a is the total zero point energy and 
Co the effective frequency given below. (•)$ denotes the 
expectation value with respect to the state |$). iJ® ff and 
-ffgp are the energy and electron-phonon coupling for the 
effective multiplet: 

5 2 
Hf = 2+ N r> where N P= Yj a k a ^ (8) 

fe=-2 

H% - -9\j\ £ Rl k nln ,clc n , s q k , (9) 

k n n' s 

where N p is the phonon number operator for the effective 
mode and 

8 8 

a 2 = ^2gl, u a = —, L) = ^uj a u 2 a . (10) 

a=l 9 a=l 

Hence, the ground state energy for the effective single 
mode model with frequency Co and coupling constant g 



Mode 


u> a (cm 1 ) 


iu a (meV) 


X a /N(0) (meV) 


3a 


H s (8) 


1575 


195.3 


22 


.368 


H g (7) 


1426 


176.8 


20 


.368 


HM 


1248 


154.7 
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H g (5) 


1099 


136.3 


12 


.325 
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708.5 


87.8 


12 


.405 


H g {2) 


430.5 


53.4 


38 


.924 


H s (l) 


270 


33.5 


21 


.868 


Eff 


581 


72.1 




1.532 



TABLE II: Frequencies and coupling constants for the vibra- 
tional modes in Cgo as taken from Manini^. The set of pa- 
rameters ori gina tes from Gunnarssoni^. The parameter sets 
used in Ref. |22| and Qjl differ marginally in the frequencies. 
The coupling strength g a and the electron-phonon coupling 
X a /N(0) are related by g\ = (6/5)A c< /(w c «iV(0)). The last 
line are the resulting parameters for the effective mode ap- 
proximation as given by equation QQ) . 



yields a variational estimate for the ground state energy 
of the multi-mode problem. Parameters cu a and g a are 
given in Tab.HTland were taken from Manini 2 -. They go 
back to photoemission measurements on gas-phase C^~ 
by Gunnarsson et a h 19 ' 23 . These parameters lead to u> = 
72.1 meV and g = 1.532. 

In the present work, we will always use the effective 
mode approximation. Hence, we will work with the 
Hamiltonian 

H c ff — Hp S + + —Hj, (11) 

where we dropped Hjj which only contributes a constant 
for fixed charge. For convenience, all energies are ex- 
pressed in terms of Co. In order to shorten the notation, 
we will omit the superscript "eff" for H p and H ep in the 
following. The effective Hamiltonian has still SO(3) and 
particle-hole symmetry. The use of the EMA is justified 
for two reasons: First, it was shown by O'Brien that it 
is a good approximation for the ground state energy and 
that multimode corrections are smalU 3 -. This will also 
be confirmed in this work when comparing the present 
results to the literature in section IVIII Second, we are 
interested in the low energy excitations of the ion. 
Generally, there are two types of excitations: rotator 
excitations involving a collective distortion correspond- 
ing to the effective mode and vibrational excitations in- 
volving "individual" modes out of the phonon spectrum. 
While the EMA is well suited for rotator excitations, it 
obviously doesn't capture vibrational excitations. While 
the former possess the energy scale Cj/(3g 2 ) « 10 meV 
(see below), the vibrational excitations lie in the range 
iv a « 30 — 200 meV. Consequently, the low-energy ex- 
citations are rotator excitations and are well described 
within the EMA. 

In what follows, it will be convenient to express the 
electron-phonon coupling term H ep = g{A\ + A ep )/^/2 
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FIG. 1: Distortion of a sphere with radial displacement 
Ar (n) oc cos a Y20 (fi) +sin a Y22 (fi) for shape parameters a = 
0, tv, tv/2. The distortion Ar has symmetry z 2 — (x 2 + y 2 )/2 
for a — 0, 7r and \fH(x 2 — y 2 )/2 for a = iv/2. The distorted 
spheres for a = 0, 7r are axially symmetric with respect to 
z-axis. The case a = it/2 is not axially symmetric, but has 
twofold axes given by the coordinate axes. 



in terms of the operator A ep defined by 



1 dfc. 



(12) 



The operators A ep and A£ p can be understood as anni- 
hilation and creation operators because they annihilate 
or create a phonon. However, they don't obey simple 
commutation relations as a k and Ofc. For this reason, the 
p®d Jahn- Teller problem doesn't have a simple analytic 
solution such as a single displaced harmonic oscillator. 



III. COHERENT STATES 

The Jahn- Teller problem under consideration involves 
distortions of the molecule which we describe by coherent 
phonon states of the type 



|q) = exp 



2 \ 

: Y qk P k 

fc=-2 / 



|o>, 



(13) 
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where |0) denotes the vacuum state and pk — 
afe)/v2 the phonon momentum operator. The vector 
q = (q~2, ■ ■ • , <&) gives the displacements of the oscil- 
lators and parameterizes the state. Generally, any 5- 
vector q can be parameterized by the overall magnitude 
q = |q|, a shape parameter a and three Eulcr angles 
= 0,7)2: 



q = qD (0) (0, 0, cosa, 0, sin a), 



(14) 



where D 2 (Q) denotes a 5 x 5 matrix with elements given 
by the Wigner D-functions D% k , (0). The effect of the 
shape parameter a is illustrated in Fig. [TJ The state |q) 



has the properties: 

[/(0)|q) 

(q|?fe|q) 
<q'|q> 
(q|# ep |q)p 



|£ 2 (©) q ), 

Qk 1 v 

71 |q> ' 

Qk, 

e -||q-qT 

A„„/(q) cl s c ri 



(15) 



The first property involves the rotation operator [/(©) 
and defines how the coherent state transforms under ro- 
tations. In the last equation, the expectation value of the 
electron-phonon coupling is determined with respect to 
the phonon state |q) which leaves an electron operator 
with matrix elements A„ n ' (q) given by 

q + \/3q 2 -V3<7_i -V3q- 2 
-V3<?-2 -a/391 qo - \/3<?2 



A(c) = f 



In order to study the eigenvalues of A(q), we invoke the 
SO (3) symmetry of the Hamiltonian. Due to this symme- 
try, the eigenvalues of A(q) and A(q') have to be equal if 
|q) and |q') = f7(0)|q) are related by a rotation. There- 
fore, in view of the parameterization (|14j) . the eigenval- 
ues of A(q) only depend on the magnitude q and the 
shape parameter a. Choosing q = (0, 0, q cos a, 0, q sin a) 
makes A(q) diagonal and yields the eigenvalues 



(A y , A z , A x ) = gq (cos 



cosa, cos 



a 



(17) 

As discussed below, the expectation value of the electron- 
phonon coupling term H ep , i.e. the eigenvalues Aj,, A z , 
A x , determine the leading term of the Jahn- Teller energy 
gain. 

In contrast, the low-energy excitations emerging from 
the rotator physics have a much smaller energy scale 
which is given by the moment of inertia of state |q) . Note 
that the moment of inertia in question is the one carried 
by the H g phonons and can be thought of as the mo- 
ment of inertia carried by a tidal wave. It should not 
be confused with the moment of inertia for overall rota- 
tions of the molecule which does not enter the present 
problem. The moment of inertia can be obtained using 
a semi-classical approac h 1 1 i 13 . It can be shown that the 
semi-classical equations of motion for the five phonon 
coordinates coupled to the ionic charge take the form 
of a quantum rotator in the limit g — > 00. This yields 
the moments of inertia. In the case q = where no 
phonons are excited, one has L|q) = and therefore no 
moment of inertia. For the state |q) acquires mo- 
ments of inertia which, due to the overall SO (3) symme- 
try of the Hamiltonian, depend only on the magnitude 
q and the shape parameter a of the distortion q. For 
q = (0,0, q cos a, 0, q sin a) one has^i: 

7T" 

6- 



, sin 2 a. cos 2 



a 



(18) 
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IV. 



The Cg ion is the simplest case among the ions. 
In particular, the Hund's rule coupling term Hj defined 
in (|U) is strictly zero and one only has to deal with the 
Jahn- Teller effect. In the following we attempt to con- 
struct variational wavefunctions starting from coherent 
states defined in (fl~3|) . The idea is to choose the electron 
wavefunction and the distortion such as to minimize the 
electron-phonon coupling H ep . This yields a state which 
gives the leading term of the Jahn- Teller energy gain. In 
a second step, the state is projected onto angular mo- 
mentum subspaces in order to investigate the rotator ex- 
citations. Finally, we improve the variational estimate 
by enlarging the Hilbert space of the trial function. 

Given the eigenvalues (fTT|) . the electron-phonon cou- 
pling is minimized by putting the electron in the ^-orbital 
(n = 0) and choosing a = 0. In this case the electronic 
levels split into a singlet with energy —gq and a dou- 
blet with energy gq/2 (neglecting spin degeneracy). This 
yields the trial wavefunction 



l*o(«)) 



(19) 



where eo = (0,0,1,0,0) is a unit vector in the 5- 
dimensional phonon coordinate space. This state is nor- 
malized and an eigenstate of A ep with eigenvalue —q/\/2. 
The expectation value of H is 



(Hi) ~2%). This is surprising as one would expect the 
Lanzcos choice to be optimal. The explanation is that 
energies are minimized with respect to q and therefore 
the procedure is not a Lanzcos expansion in the proper 
sense. From these findings we deduce the following rule: 
Given a trial function |1), we achieve a good improve- 
ment by adding state |2) = A| p |l) to the Hilbert space 
of trial functions. Note that A\ p is a scalar and there- 
fore |1) and 1 2) have the same symmetries. In addition, 
using A\ p instead H to create a second state yields sim- 
pler wavefunctions. Below, we make extensive use of this 
rule. 

The state considered above is not an angular mo- 
mentum state, but it can be understood as a rotator at 
rest. Its moments of inertia are given in (|18p . For a = 
and q = g we find / — I x = I y = ig 2 and I z = 0. There- 
fore we expect a rotator spectrum given by L(L+1) / (6g 2 ) 
which is on a smaller energy scale than the Jahn- Teller 
energy gain of the order g 2 . The rotator takes up 2 de- 
grees of freedom of the 5-dimensional phonon space. The 
other three degrees of freedom are vibrations of the ro- 
tator. This picture emerges also when treating the p ® d 
Jahn- Teller problem semiclassically as done in Ref. Ulllig . 
These references show that the asymptotic rotator spec- 
trum in the g — * oo limit becomes 



E 



9_ 

2 



L(L + 1) 
6g 2 



(23) 



(20) 



which is minimal for q — g and which yields an upper 
bound for the ground state energy: 



E 



N=l 



2 



(21) 



The energy —g 2 /2 is the leading term of the Jahn- Teller 
energy gain of the Cg ion. 

The variational estimate of the ground state energy 
above can be improved by enlarging the Hilbert space 
for the variational wavefunction. We consider three dif- 
ferent Hilbert spaces spanned by the following choices of 
wavefunctions 



(i) l*o) 
(**) l*o) 
(Hi) |* ) 



H|*o) 
4pl*o> 



(22) 



4pl*o> 



Choice (i) corresponds to a Lanzcos step. In choice (ii), 
A\ p creates a phonon excitation. Choice (Hi) allows for 
both. The basis spanned by each choice depends on the 
variational parameter q and the variational ground state 
energy is obtained upon minimization of the lowest eigen- 
value with respect to q. Obviously, choice (Hi) must 
yield the lowest estimate. It turns out that choice (ii) 
and (Hi) yield almost the same energies (rel. diff <0.2%) 
whereas choice (i) is somewhat higher (rel. diff to choice 



The first term is the leading Jahn- Teller energy gain. The 
second term is the zero-point energy of the three remain- 
ing decoupled oscillators and the last term are the rotator 
excitations. 

Before going into the details of the projection tech- 
nique, we need to investigate the symmetries of \^o) in- 
order to know which projections are non-zero. First we 
note that L z \^q) = because |*o) is constructed from 
operators cj» and oj for which the z-component of angu- 
lar momentum vanishes. Hence, only M — projections 
are non-zero. In addition, due to the electron operator 



|\&o) is odd under a 7r-rotation around the y-axis. 
Therefore only projections with odd total angular mo- 
mentum L are allowed as can formally be shown using 
the second projection operator property given in (|49[) . 
These findings agree with the literatu re 11 ' 13 . 

Projection operators onto angular momentum sub- 
spaces are given by 



Qmk 



2L + 1 
8ir 2 



d@D L MK (Q)U(Q), 



(24) 



where U (0) is the rotation operator and 
Dj^ K (@) the real Wigner D-function. J dQ = 

/ Q 27r d(j> Jq d9 sin 9 J* dj denotes the integration over the 
Euler angles (<f>, 9, 7) . The prefactor serves for proper 
normalization and arises from the orthogonality relation 
of the Wigner-D functions^. These projection operators 
have the following properties with respect to normalized 
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angular momentum states \LM) 

(LiM l \Q% TK \L 2 M 2 ) = S LLl S LL2 8 MMl Skm 2 - (25) 

The lowest state is expected for L = 1 which leads to the 
trial function 

\Po(q)) =Qool*o(9)) =Q 1 00 cl ] \qe ). (26) 

We have adopted the letter P, as in atomic physics, to 
indicate that the state has total angular momentum L — 
1. The spin degeneracy is not indicated because it is 
always 2 in the case of one electron. Note that this trial 
function was already proposed in Ref. \vA but without 
having q as a variational parameter. 

Since |\I>o) 1S an eigenvector of A ep , \P ) is an eigen- 
vector as well because A ep is a scalar and [A ep , Qqq] = 0. 
Hence, we have (Po\H ep \P ) / (P \Pq) — ~gq- The calcu- 
lation of the expectation value of the phonon number N p 
is more involved. The norm of the wavefunction \Pq) is 
given by 

(P |P ) = (*o|Qj l*o) (27) 

= JL J d&D l 00 (e)(* \u(e)\* ), 

In the first equation we used the projection operator 
property (<2o ) 2 = Qho- As discussed in the introduc- 
tion, this step is crucial because eliminating a projection 
operator eliminates an integration over Euler angles. In 
the second equation of (12~7| the definition of Qq is sub- 
stituted. The integral above involves the matrix element 



(*o 1*7(6) |*o) = cos(9exp 



^2-2 

- - o sin i 
4 



(28) 



which can be calculated using the rotation rule for tensor 
operators (|A8|) as well as properties (|T5j). The expecta- 
tion value (Po|iVj,|Po) is evaluated likewise and involves 
the matrix element 



(*o\N p U{e)\* ) 



(29) 



1 

— cos( 
2 



- cos 

2 2 



exp 



^9-2 

- - q sin i 
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FIG. 2: Upper panel: Ground state (L = 1) energy and ro- 
tator excitations (L = 3, 5, 7) for C§ as a function of the 
coupling strength g and relative to the asymptotic energy 
—g 2 /2. The dashed lines are the variational energies for 
states as given by (|26p . but with an appropriate projector 
Qoo- The solid lines are the energies of the improved vari- 
ational approach. The dotted lines are asymptotic behavior 
given in (|23p . Lower panel: Energy of the ground state of 
C^q. Solid lines (a,b): Present results, the lower curve (b) 
is the improved version. Long dash (c): Numerical result of 
O'Brien''. Short dash (d): Variational wavefunction of Dunn 
et al. (curve (b) in Ref. \lw . Dots (e): Approximate analyti- 
cal result of Chancey^. Dash-dots (f): Asymptotic behavior 



given in (123 



where it should be noted that [N P ,U(Q)] = because 
N p is a scalar. The matrix elements (|28|) and (|29p don't 
depend on the Euler angles <j) and 7 which leaves one 
integration over 9. Substituting t = cos 9 we find 



(Po\N p \Py 

(Po\Po) 



[¥ 2 



-!? 2 (i-t 2 ) 2 {q> ' 

(30) 

Note that the integrals can be expressed in terms of er- 
ror functions. The last equality defines the function h(q) 
which varies smoothly from h(0) = 2/5 to h(oo) = 1. 
Putting the different parts together, the expectation 
value of H becomes 



E N=1 (g,q) 



Pn\H\P 







(P \Po) 



gq. (31) 



The only difference to the energy found above is the 
factor h(q) which renormalizes the phonon energy q 2 /2. 
This additional energy gain is due to the "delocaliza- 
tion" of the distortion in the projected state. Minimizing 
Eo(g, q) with respect to q for a given coupling strength g 
yields the variational ground state energy shown in Fig. [2] 

Higher rotator excitations are constructed by replac- 
ing the projector Qq in (|26[) with a projector Qq on a 
higher angular momentum space. This results in replac- 
ing the Wigner D-function Dq Q in integral (|27|) with Dq Q . 
Note that the Wigner D-functions £>oo( ) = Pl(cos9) 
are given by the Legendre polynomials Pr,. Results for 
L = 3, 5, 7 are shown in Fig. (5) As can be seen, the 
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asymptotic behavior for g — > oo differs substantially from 
the semiclassical result ([23"]). In order to improve these 
energies we follow the rule suggested above. That is, 
we use a 2-dimensional Hilbert space spanned by \Pq) 
and ^4| p |Po) and minimize 25 the lower eigenvalue with 
respect q. Higher angular momentum excitations are 
treated likewise. Note that At is a scalar and there- 
fore [Qm K ,A\ ] = 0. Results are shown in Fig. O The 
improvement is substantial and the asymptotic behavior 
for g — ► oo fits well the semi-classical result (|23l) . 

Comparison of the present result for the ground state 
energy of Cg with other works are shown in the lower 
panel of Fig. [2] Our improved result fits well the numeri- 
cal calculation of O'Brien^ in the low coupling regime. 
It deviates for higher couplings where the energies of 
O'Brien are somewhat lower. This might be due to the 
fact that our approach is strictly variational whereas in 
Ref. [9| the matrix elements for states with a large num- 
ber of phonons are extrapolated which may lead to non- 
variational energies. In the limit g — > oo, our result 
approaches the asymptotic behavior (|23p smoothly from 
above. 



A new aspect of Cqq is the non-trivial Hund's rule 
coupling Hj given in Q . The main features of the com- 
petition between Jahn- Teller effect and Hund's rule cou- 
pling can be observed on the level of unprojected states 
which we discuss first. As before we start with the follow- 
ing state which minimizes the electron-phonon coupling 
H, 



ep- 



l x *o(g)) 



c OT c o|l<? e o> 



(32) 



The upper index indicates that the state is a spin singlet 
with spin degeneracy 1. This state is normalized and an 
eigenvector of A ep with eigenvalue \J~2q which is twice 
bigger than for Cg due to presence of 2 electrons. The 
expectation value oiHj with respect to l^o) is 2J which 
implies, in view of the multiplet energies given in Tab. [H 
that I 1 'to) is not an eigenstate of Hj. In fact, Hj couples 
I ^o) t° another state 



r*i(<7)) 



V2 



L — IT — 14. 



ke ), 



which also has spin 5 = and total angular momentum 
Lz^^i) = 0. It is an eigenstate of A ep with positive 
eigenvalue q/V%- The two states l^o) and l^i) form 
a basis in which A ep is diagonal but Hj is not. Eigen- 
states of Hj are obtained by the following orthogonal 
transformation: 

1*5) = ^ = -^£ 4Akeo>, 

v v n— — 1 



I 1 * 



D 



-v / 2| 1 ^o) + | 1 ^i) 
V3 



(34) 



The corresponding eigenvalues are 4 J and J respectively, 
as given in Tab. HI These two states are denoted with 
the lower indices 5 and D because they are eigenvectors 
of the total electron angular momenta L 2 with angular 
momentum L e — and L e = 2 respectively. The Hamil- 
tonian in the basis ( 1 4's, *&d) takes the form 



5 <? f U/u) V2gq 
2 2 V ^91 -91 + J /u 



(35) 



Minimizing with respect to q yields the ground state en- 
ergy 



E, 



N=2 



J 



H + - - g*f — 

1 1 iog z 



J 



(36) 



where the function / decreases monotonically from 
/(0) = 2 to /(oo) = 1/2. Hence, in the absence of 
Hund's rule coupling (J = 0), the Jahn- Teller energy 
gain is 2g 2 . This is 4 times bigger than in Cg because 
the electron-phonon coupling is doubled due to the pres- 
ence of two electrons. In the case of dominating Hund's 
rule coupling, i.e. J 3> wg 2 , the Jahn- Teller energy gain 
is reduced to g 2 /2 but not entirely suppressed because 
^^olHep^^D) = -gq. Note that this finding differs 
from the general belief that strong Hund's rule coupling 
completely suppresses the Jahn- Teller effect. 

If we don't restrict our view on the 5 = sector, then, 
of course, the spin triplet state will be favored for large 
enough and positive J (see Tab.|T|. Hence, there is a level 
crossing between low and high spin state which depends 
on the parameters g and J. The Jahn- Teller problem in 
the 5=1 sector is equivalent to the Jahn- Teller prob- 
lem of Cg~ due to particle-hole symmetry. This is obvi- 
ous when looking at the case of maximal spin S z = 1, 
where the spin-up states are occupied by two electrons 
and one hole. Therefore, the energy of the triplet state 
on this level of approximation is E = 5/2 — ,g 2 /2 — J/ui, 
This energy can be compared to the energy of the sin- 
glet state (1551) . One finds that the level crossing occurs 
at J/(uog 2 ) = 0.5284. This criterion is only little modi- 
fied when going to projected variational states below (see 
Fig. E]) . 

As before, we calculate rotator excitations using pro- 
jection operators. The phonon coherent state in the 
(33) present case is the same as for Cg . Therefore we ex- 
pect the same rotator physics in the strong coupling limit 
g — > 00. The moment of inertia of the rotator is / = 3<7 2 . 
In the absence of Hund's rule coupling, where q = 2g, we 
obtain / = 12g 2 which implies that the energy scale of the 
rotator excitations in Cg^ is four times smaller than in 



and J = 



C 60 . The full asymptotic behavior for g - 
is again obtained from the semiclassical approach 13 : 



E 



n=2 



~2. 9 2 



1 L(L + l) 



12 5 2 



24.g 2 



(37) 



As in C 60 , the unprojected states are annihilated by L z 
which means that only M — projections are allowed. 
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On the other hand, the states are invariant under ir- 
rotations around the y-axis which requires L to be even 
(see properties (j49j) of the projection operator). 

In order to investigate the rotator physics in the pres- 
ence of Hund's rule coupling we use the following pro- 
jected states: 



\ 1 X 8 ) = Q, 



L il 



\ 1 X D )=Q 



L il 



D 



(38) 



where X = S,D,G,I, . . . stands for the letter denot- 
ing the total angular momenta L = 0, 2,4, . . . The two 
states l^s) and I^d) are eigenstates of Hj with differ- 
ent eigenvalues and so are the two states defined in l|38p 
because [Hj, Qq ] = 0. Therefore, they form an orthogo- 
nal basis. The calculation of the various matrix elements 
follows the procedure described above. Within basis ([38]) 
and for a given L, the Hamiltonian has the matrix ele- 
ments 



('Xs^Xs) 
CXnlH^X 



D 



I^Xd^Xd) 
^^X S \^X S )^X D \^X D ) 



where the function (q) is defined in terms of Legendre 
Polynomials Pj,: 




2L+1 



dtP L (t) [P 2 (t)] n e 



(40) 



The energies which result from minimizing the lower 
eigenvalue are shown in Fig. [3] for J = and J = ui. 
As expected, J leads to an inversion of the L = and 
L = 2 levels for small enough g. For J — Q and 5 = the 
L = level has energy 4.5<D. This state corresponds to 
an electronic -D-state with one phonon excitation. The 
pure electronic S'-state with no phonon excitation has 
higher energy 6.5<D. For large enough g the two spec- 
tra become very similar except for an overall energy shift 
~ 2 J. This corresponds to the expectation of H j with re- 
spect to state | ^o) which minimizes the electron-phonon 
coupling. 

As can be seen from Fig. [3l the asymptotic energies 
as calculated above overestimate the rotator excitations 
substantially. Therefore we improve the variational ap- 
proach as discussed above. That is, the 2-dimcnsional 
Hilbcrt spanned by (f3"5|) is enlarged to a 4-dimensional 
space using, in addition, the two states Al^Xg) and 
A\ ^Xd). The improved energies^ are also shown in 
Fig. [3J The asymptotic behavior agrees well with the 
semiclassical result (|37|) . 

For big enough Hund's rule coupling J, a level crossing 
occurs such that the S = 1 spin triplet state becomes 
the ground state. The crossing between the spin singlet 
and spin triplet ground state defines a line in the (g, J) 




FIG. 3: Upper panel: Rotator states L — 0, 2, 4, 6 of C 6 q in 
the singlet subspace (S = 0) and for J = 0. The levels are 
plotted as a function of the coupling strength g and relative 
to the asymptotic energy —2g 2 . Dashed lines are the results 
using basis (|38p . Solid lines correspond to the improved ap- 
proach and the dotted lines are the asymptotic behavior given 
in (|37|) . Lower panel: As in the upper panel but with J = Q. 



parameter space as shown in Fig. [5l Note that singlet- 
triplet crossing occurs always for smaller J than the L = 
to L = 2 crossing in the spin singlet sector. The energy 
of the triplet state is obtained by subtracting J from 
the ground state energy of C^q. The crossing line was 
calculated using improved variational approaches. As can 
be seen in Fig.0 the criterion J/(ujg 2 ) = 0.5284 derived 
above for the unprojected states becomes correct in the 
large g limit. For g — * the line ends at J/(ujg 2 ) = 3/4 
which can be shown using the perturbative results^ for 
small g. 



VI. 



In order to find a state which minimizes H ep we 
start with the distortion q = q(0, 0, cos a, 0, sin a) and 
put two electrons in the x-orbital and one in the z- 
orbital. According to (fl~7|) , the expectation value of H ep 
is 2Aj; + A z = —^/igqsa\a which is minimal for a = tt/2. 



9 



This yields the spin-1/2 state 
| 2 *o(<z)> 



(41) 



where the upper index denotes the spin degeneracy and 
e 2 = (0,0,0,0,1). State \ 2i S?o) is an eigenstate of A ep 
with eigenvalue — ^3/2q. However, it is not an eigenstate 
of Hj. Hj couples | 2 ^o) to another state 



*i(?)> = cL c -iT c -nl ( ? e 2>> 



(42) 



which is an eigenstate of A ep with positive eigenvalue 
\/3/2g. Eigenstates of Hj are given by 



2 *o) 



V2 



I 2 *, 



V2 



(43) 



These states are also eigenstates of the total electron an- 
gular momentum as indicated by the lower indices P and 
D. Within basis (14"3"1), the Hamiltonian has the form 



( 2J/Cu V3gq \ 
2 + I V3gq J ■ 



(44) 



Minimizing the lower eigenvalue with respect to q yields 



EST 6 = 



2 



J 



J 



6u 2 g 2 



J > iujg 2 



< J < 3ug 2 



(45) 

For J = the Jahn- Teller energy gain is —3g 2 /2 and 



the 



somewhat reduced compared to C 6 q . For J > 3Qg 
Jahn- Teller effect is completely suppressed. 

As in Cgp , there is a high spin state (S = 3/2) 
which is favored by the Hund's rule coupling. In this 
state, the electrons have parallel spin. Hence, each p- 
orbital is occupied by one electron and no Jahn- Teller 
coupling is possible. According to Tab. [H the ground 
state energy of the S = 3/2 state is simply given by 
5/2 — 3 J. Using energy (|45[) one finds the criterion 
J/(cDg 2 ) = 3(4- yT5) = 0.381 for the low-spin/high-spin 
crossing. 

In order to proceed with the projection, we first ana- 
lyze the moments of inertia. According to relation (fTS]) , 
there are two different moments of inertia I\ — I x — 
Iy = q 2 and h = I z = 4g 2 . The fact that I x = I y and 
I z ^ suggests that the rotational degrees of freedom be- 
have similar to a symmetric top 11 ^ 13 . The solution of the 
quantum mechanical symmetric top is well known (for an 
early review see Ref . \2w . The symmetric top has three 
rotational degrees of freedom, the Euler angles, and three 
conserved quantities: the rotational frequency around 
the principal axis -Z3 of the top (quantum number K ) and 
the angular momentum corresponding to the precession 
of this axis (quantum numbers L,M). The two remain- 
ing degrees of freedom become decoupled oscillators, i.e. 
motions of the top axis. The spectrum of the symmetric 




FIG. 



4: Upper panel: 
,3- 



Rotator states (L K) = (10), (2 - 
2), (4 - 4), (3 2) of Cgo in the S = 1/2 sector and for J = 0. 
The levels are plotted as a function of coupling strength g 
and relative to the asymptotic energy — 3<? 2 /2. Dashed lines 
are the result using basis (|5ip . Solid lines correspond to the 
improved approach and the dotted lines are the asymptotic 
behavior given in (|46[) . Lower panel: As in the upper panel 
but with J = Q. 



top26 i s given by (If 1 L(L+1) + [I3-J1 1 }K 2 )/2. This 
spectrum is recovered for the low energy excitations of 
Cqq in the limit g — > 00 where q = \/3g, I\ = Sg 2 and 
7-5 = 12g 2 . Within the semiclassical approach 13 one finds 
the asymptotic behavior 



E^ A = - 



2 



24 5 < 



+ 



L{L+l)-lK 2 
6g 2 



(46) 



Eigenfunctions of the top are given by the Wigner-D 
functions D^ IK . Since we chose the distortion q such 
that the 7 3 -axis corresponds to the z-axis, states with 
good quantum numbers L, M, K are given by the projec- 
tion Qmk- 

As above, there are selection rules given by the sym- 
metry. In contrast to the previous cases, the states \ 2 ^) 
defined above are not axially symmetric and therefore 
projections with K ^ are possible. The remaining 
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symmetries of these states are 

t/(0,0^)| 2 *) = | 2 *), 

t/(0,7r,0)| 2 vK) = -I 2 *), 

SC/(0, 0,7r/2) I 2 *) = | 2 *). 



The matrix elements are given by: 



(47) 



( 2 X« \H\ 2 X*} 

( 2 Xf | 2 X£) 
( 2 X«\H\ 2 X«) 



where the arguments of the rotation operator U(<f>,9,"f) 
are the three Euler angles. The last symmetry involves 
the particle-hole transformation S. This symmetry is 
special to CgQ where the number of holes is equal to the 
number of electrons. The operator S is defined by 



( 2 Xg\H\ 2 X. 



5 

2 + 
5 

2 + 

x/3 
2 



q 2 G 



LK 



2 Nl K 



2J 

LU 



q 2 G 



LK 



(52) 



gq 



N D 

LK 



p 

LK 



"Ofe, 



Sc ns S 



(48) 



The definition implies that S commutes with the angular 
momentum, i.e. [L, S] — 0, and therefore [U(Q),S] = 
[Qmk> = 0. Selection rules for the quantum numbers 



LK can be derived using the symmetries (I47|) and the 
following properties of the projection operators 



where 
NE K (q) 

nEk(i) 
GlAi) 



G D LK {q) 



21 + 1 

2/ + 1 
"S^r 2 " 
2Z+1 

8tt 2 
2Z + 1 

8tt 2 



dQD L KK Dl e- q ^ l - D ^\ 
dQD L KK D*_ 2 _ 2 e-^- D ^\ (53) 
dQD KK Dl D 2 2 e-^ l - D ^\ 
deD KK D 2 _ 2 _ 2 D 2 2 e-^^- D ^. 



Q^ k U(0,0,k) = 
Qi IK U(0,n,0) = 

Q L MK SU (0,0^/2) = 



K Q L 



(-1 

*k(-1) L 



)mk> 



(49) 



(-1)^ SQi IK for K even. 



The first property implies that K is even. The second 
property implies that K > for odd L and K < for 
even L due to the definition (|A3[) of . Note that K = 
states are forbidden for even L. These rules agree with 
the literature^. The last property, applied in the second 
equality below, yields the rule for the orthogonality of 
projections: 



In the expressions above, the integration over the Euler 
angles </>, 7 is not trivial. As can be seen from (|A7|) . 
the Wigner D-functions D^ K depend on 0,7 through 
cos[A'((/j±7)]. Therefore, the integration over </j, 7 can be 
carried out in terms of modified Bessel functions /„ using 
the new integration variables /i = 4> + 7 and v = <p — 7. 
For example, N^q becomes 



dtt 1 



q 2 (l + t) 2 



8 



Io 



<7 2 (l-i) 5 



(.54) 



2 ^\Qk,m,Q 



Li 



K 2 M 2 ^M 1 K 1 



(-l)-^J 



5LiL 2 °J\/i M 2 



[ 2 MQ l k\ Ki \ 2 ^) 



The remaining integrals are numerically evaluated. The 
lower eigenvalue is then minimized with respect to q in 
order to find the variational ground state energy. Fig. 0] 
shows the resulting energies for J = and J = Q. The 
asymptotic spectrum (|4"6")l is also shown in Fig. 2] and dif- 
fers from the energies of the present calculation. Again, 



L 1 L 2 °M 1 M 2 



^\QfC k | 2 ^)- (50) t ^ s l& i m P rovec i by adding two more states, A\ p \ 2 X 



K\ 



This implies that projections with equal L and M are 
only orthogonal if (K\ + K 2 )/2 is odd. With these se- 
lection rules and the spectrum of the top given in (|46|) . 
we find that the ground state and lowest excitations have 
quantum numbers (LK) = (10), (2-2), (4-4), (32), (30). 
The next higher state in energy is (6—6) which, according 
to (|50p. is allowed to mix with (6 —2). 

In order to calculate rotator excitations, the basis 
states (|43| are projected: 



\ 2 X*) 



i 0K I 



\ 2 x%) 



Q 



OK 



(51) 



Again, X denotes the total angular momentum and K is 
the quantum number for the rotation around the princi- 
pal axis of the top. These two states are eigenstates of 
Hj with eigenvalues 2J and respectively. The calcu- 
lation of the expectation value of H proceeds as before. 



and A\ \ 2 Xp), to the Hilbert space (see Fig. [4 

As discussed above, a level crossing from a low-spin 
to a high-spin state occurs in Cg^~ for large enough J. 
The corresponding line in the (g, J) parameter space is 
shown in Fig. O and was calculated using the improved 
variational approach. The criterion J/(ujg 2 ) = 0.381 de- 
rived above becomes correct in the large g limit. For 
g — ► 0, the line ends at J /(tug 2 ) = 3/4 which can be 
shown using perturbative results^ for small g. 



VII. RESULTS FOR C 60 

In this section we calculate ground state and excita- 
tion energies for parameters specific to Cgo- As shown 
in section [TH the ground state energy in the effective 
mode approximation (EMA) is a variational estimate for 
the ground state energy of the full multi-mode model. 
Ground state energies and lowest excitations obtained 
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N 



E(P) 



E(F)-E(P) 
E{H) - E{P) 
E(L) - E{P) 



E( L D)-E( l S) 
ECG)-ECS) 
£( 1 J)-S( 1 5) 



E( a P)-ECS) 
W 



E( 2 P U ) ~ 
E( 2 D 2 )~E( 2 P°) 
E{ 2 G l )-E{ 2 P°) 
E( 2 F 2 )-E{ 2 P°) 



E(' l S)-E( 2 P u 



J = 
Ref. 22 Present 



-139.6 -132.8 



26.8 
53.9 



40.9 
87.4 
137.3 



-126.7 -141.2 



-405.9 -406.8 



13.2 
38.4 
74.8 



29.5 
56.0 



266.3 273.9 



329.3 345.1 



-342.9 -335.7 



5.0 
28.1 
35.4 



342.9 335.7 



-126.0 -142.2 



<:), 



j 



-132.8 -137.6 



40.9 
87.4 
137.3 



-20.6 
-37.4 
-52.0 



-5.6 



-167 1.75 



-271.2 -442.2 1.75 



9.6 
29.4 
56.2 



-9.8 
-31.4 
-55.6 



0.003 
0.010 
0.018 



66.8 304.6 -2.75 



141.6 386.5 -2.63 



-268.0 -360.3 0.87 



4.8 
28.0 
35.6 



-4.2 
-22.6 
-24.8 



-0.002 
0.002 
0.007 



51.7 360.3 -3.87 



-6.4 -163.; 



1.76 



TABLE III: Ground state energies and low-energy excitations 
for C^) - (N = 1, 2, 3) as calculated by the effective mode ap- 
proximation (EMA). [/'(AO = E (N~1)+E (N+1)-2E (N) 
is the contribution of the Jahn- Teller and Hund's rule coupling 
Hamiltonian H' — H p + H ep + Hj to the effective on-site re- 
pulsion 19 U. Improved variational states with g = 1.532, 
Gj = 72.1 meV, J = (columns 3,4) and J = u) (columns 
5-7) are used. All energies are in meV. The zero-point energy 
ojq = ~}2 a (5/2)ujct is neglected. Column 3 is the exact diago- 
nalization result from Ref|22j. Columns 6,7 are the derivatives 
with respect g and J. 



from the EMA using improved variational approaches are 
given in Tab. IIII1 Concerning the phonon-related param- 
eters g and Q, there is a consensus in the literature to use 
the parameter set given in Tab. HT1 which originates from 
photoemission experiments on C^~ in gas phased. This 
parameter set yields g = 1.532 and uj = 72.1 meV for the 
EMA. On the other hand, there is much less consensus 
concerning the Hund's rule coupling J. Below, we will 
first discuss the case J = where exact diagonalization 
results are available^. Subsequently, we will determine 
the parameter J > through the singlet-triplet gap in 
Cqq which can be measured experimentally. 

We first consider the case J — 0. As can be seen from 
Tab. IIII1 the relative errors between the present ground 
state energies and the exact diagonalization results^ are 
4.9%, -0.2%, 2.1% for N = 1,2,3 respectively. Hence, 
the present results agree well which confirms the valid- 
ity of the EMA. The agreement is best for N — 2 where 
the variational energy is in fact below the exact diago- 
nalization result. This is due to the truncation of the 
phonon Hilbert space in the exact diagonalization ap- 
proach which makes it variational as well. Hence, exact 
diagonalization tends to overestimate the ground state 
energy for a large number of excited phonons which is 
the case for N = 2 where the Jahn- Teller energy gain is 
largest. The present approach based on coherent phonon 




0.4 Low Spin 

0.3 ' ■ — 

0.5 1 



® Bulk 



1.5 

g 



0.381 

2T 



FIG. 5: Boundaries between the low-spin and high-spin sec- 
tors in the (g, J) parameter space. The upper line is for N — 2 
and indicates the 1 S- 3 P level crossing whereas the lower line 
is for N = 3 indicating the 2 P°- 4 S level crossing. Both lines 
end at J/(u>g 2 ) = 3/4 for g — *• 0. For g — ► 00, the up- 
per line approaches the constant value 0.5284 whereas the 
lower line approaches 3(4 — VTE) = 0.381. As discussed in 
the text, realistic parameters for bulk Ceo are g = 1.532 and 
J = Q = 72.1 meV which is indicated by the cross. 



states doesn't suffer from this truncation. The lowest 
excitations are also given in Tab. IIIIl For N = 1, the ex- 
act diagonalization yields lower energies than the present 
approach opposite to = 2. This is again due to the 
truncation effect discussed above. In addition, the exact 
diagonalization yields also the vibronic excitations which 
are not captured by the present approach. As discussed 
in Ref. [22j, the two lowest levels for N = 1 are the L = 1 , 3 
rotator states whereas the third level are L — 2 vibronic 
excitations^. In the case N = 2, the L = 0, 2, 4, 6 rota- 
tor states are lowest in ener gy, f ollowed again by L = 2 
vibronic excitations — . Ref. [22J doesn't provide excita- 
tion energies for N = 3. However, as can be seen from 
Tab. IIIIl the low-energy rotator excitation for N = 3 are 
smaller than those for AT = 2. This suggests that the 
four lowest levels of N = 3 are also pure rotator states. 

In the literature, estimations for the value of J dif- 
fer largely. Theoretical values range between J = 15 — 
300 meV (see Ref. H and references therein). Experi- 
mentally, J is not directly accessible. However, the low- 
spin/high-spin gap for A^ = 2,3 can be measured by 
various means. As can be seen in the last column of 
Tab. IIIIl this gap depends strongly on J. Below we de- 
termine the value of J using experimental values for the 
low-spin/high-spin gap. There is a consensus that iso- 
lated C$>~ (N = 2,3) ions arc in the low-spin state (see 
Ref. [I?] and references therein). However, there was a 
controversy^! on whether the low-spin/high-spin gap is 
very small (below one wavenumber— ) or rather of the 
order of 600 wavenumbers^. Recently, this problem was 
carefully reconsidered and it was shown that activated 
behaviors of C^ - which were observed so far and used 
to determine the gap are in fact due to C120O impuri- 



tie; 



,30,31 



However, the work clearly reconfirms that iso- 
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lated C 60 ions (N — 2,3) are in the low-spin state. In 
view of Fig. [5] this implies that J/ (g 2 Co) < 0.5 which 
yields the upper bound J < 85 meV for isolated C^~. 
Measurements of the low-spin/high-spin gap exist for 
Ceo bulk materials, in particular for K4C60 which is a 
non-magnetic insulator. The Cqq ion in this material 
is equivalent to the CgQ ion by particle-hole symmetry. 
Hence, the low-spin/high-spin gap should correspond to 
the singlet-triplet gap A$t of CgQ . This gap is observed 
in magnetic susceptibility 3 ^, and spin relaxatio n 33 i 34 i 35 i 36 
measurements on K^Cgo- The magnetic susceptibility 
and the spin relaxation scale with the thermal occupa- 
tion of the triplet 3 P state which shows an activated 
behavior with a gap A$p ~ 50 — 100 meV. Using the 
J-dependence of the gap as given in Tab. IIII| one de- 
duces J « 60 — 80 meV in agreement with the upper 
bound 85 meV found above. For convenience we suggest 
J = Co — 72.1 meV which leads to low-spin/high-spin 
gaps of 66.8 meV and 51.7 meV for N = 2,3 respectively. 
Energies for J = Q are given in Tab. IIIII together with 
derivatives with respect to g and J. The J-dependence 
in the low-spin sector of N = 2, 3 is non-trivial. However, 
to a good approximation, levels are shifted linearly and 
in parallel for a given N . Therefore, excitation energies 
in the low-spin sectors depend little on J as is confirmed 
in Tab. [En] 

Comparing the ground state energies for the cases 
J = and J = Q in Tab. IIIII shows that the Jahn- 
Teller effect is partly counterbalanced by the Hund's 
rule coupling. This observation is particularly relevant 
for the corrections to the on-site repulsion U. Gener- 
ally, the main contribution to the effective on-site re- 
pulsion is the isotropic Coulomb repulsion Hjj which 
we separated from the Hund's rule coupling in Hamil- 
tonian (JTJ). In addition, there is a second contribu- 
tion^ U'(N) = E (N - 1) + E (N + 1) - 2E( N) 
which is due to the different ground state energies of 
H' = Hp + H ep + Hj for different occupation numbers 
N. It was argued that this contribution is not negligible 
and may explain why compounds with average occupa- 
tion number N = 2, 4 are insulating whereas compounds 
with N = 3 are mostly metallic. Indeed, for J = we 
have U'(2) — U'(3) « 0.5 eV which is important compared 
to U = 1-2 eV. For J = Co, the difference J7'(2)-f7'(3) is 
reduced to 0. 14 eV which is an order of magnitude smaller 
than U. Hence, including J reduces U' significantly. 



VIII. CONCLUSION 

With the present approach, variational wavefunctions 
for the ground state and rotator excitations of Cq ~ 
ions are constructed semi- analytically. The Jahn- Teller 
physics, where the EMA is used, and the Hund's rule 
coupling are treated on the same level. The strict use 
of the SO(3) symmetry and projection operators allows 
for an efficient formalism. In this formalism it is evident 
that any scalar operator, such as the Hamiltonian itself, 



commutes with any projection operator. Thanks to this 
property, only one integration has to be done numeri- 
cally in the final expressions of the expectation values. 
This is a major achievement over previous approaches 
and enables the present approach to go beyond previous 
results. 

In a first step we calculate ground state energy and ro- 
tator excitations for the three distinct cases N — 1,2,3 
whereby each case has some additional complication com- 
pared to the previous one. The simplest case is N = 1 
and served to explain the projection operator technique 
in detail. We find low-energy excitation in agreement 
with previous works. The correct asymptotic behavior 
is recovered when the improved version with an enlarged 
Hilbert space is used. The main challenge for N = 2 is 
the additional Hund's rule coupling. We investigated the 
competition of Jahn- Teller effect and Hund's rule cou- 
pling on the level of both, projected and unprojected 
states. We find, somewhat in contradiction to the general 
picture, that strong Hund's coupling doesn't completely 
suppress the Jahn- Teller effect, but rather reduces the 
effective electron-phonon coupling constant by a factor 
of 2 within the 5 = sector. Of course, strong enough 
Hund's rule coupling favors the S = 1 state. We calculate 
the separation between the low-spin and high-spin sector 
in the complete {g, J) parameter space. The difficulty of 
N = 3 lies in the fact that the unprojected state minimiz- 
ing the electron-phonon coupling is not anymore axially 
symmetric. The problem therefore becomes similar to a 
symmetric top and states involve a third quantum num- 
ber K. Using the symmetries of the unprojected state 
we deduce the allowed values for the quantum numbers 
LK in agreement with previous findings. A new result 
is, that two states with L\ = L2 and odd (K% + iG)/2 
are allowed to mix. The evaluation of matrix elements for 
N = 3 is more complicated and involves Bessel functions. 
Nevertheless, only one numerical integration is required. 

Using the results of the previous sections, we calculated 
ground state energy and lowest excitations in agreement 
with exact diagonalization results. In addition, we give 
a thorough discussion of the parameters specific to Cgo- 
Whereas there is a consensus on the value of the electron- 
phonon coupling, there is much uncertainty on what con- 
cerns the Hund's rule coupling. We use the present re- 
sults to make a connection between the Hund's rule cou- 
pling constant and the low-spin/high-spin gap which is 
experimentally accessible. This allows us to pin down the 
Hund's rule coupling constant to J = 60 — 80 meV. Using 
J = Q = 72.1 meV to calculate the ground state ener- 
gies, we find that the finite Hund's rule coupling partly 
counterbalances the Jahn- Teller energy gain and that the 
ground state energies for the cases N — 2,3,4 become 
almost equal. Therefore, the contribution to the on-site 
repulsion arising from the Jahn- Teller effect is substan- 
tially reduced when including the Hund's rule coupling. 

The authors thank T.M. Rice, C. Helm and I. Milat 
for fruitful discussions. This work has been supported by 
the Swiss Nationalfonds, by the NCCR MaNEP and by 
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Symmetry: 


SO(3) 


h 


SO(3) 


h 


IR: 


L = 1 


tl u 


L = 2 


H 9 




M = —1 


It 

y 


M = — 2 


z 






M = 


z 


M = — 1 


X 






M = 1 


X 


M = 
M = 1 
M — 2 


a/37861- a 

y 

a/5786 1 + a 


/578e 
/378e 



TABLE IV: Relationship between the components of the L = 
1,2 IR's of SO(3) as used in the present work (columns 2 
and 4) and the components of the ti u and H g IR's of the 
icosahedral symmetry Iu as defined in Ref. |3?] (columns 3 and 
5). Note that the notation used in columns 3 and 5, i.e. the 
letters x,y,z,9,e denoting the different components, follows 
Ref. [H 



same way: 



E A * 



i 

E 

m— — 1 



A„ m (-1)' 



(A5) 

The operators a£ and at are given by the same rules. 

|0) = 

|0) whereas the electrons created by c\, s have y, 
-1,0,1. Using the defini- 



Note that cj, ls creates an electron such that L z c ms 

mc, 



or x symmetries for n 



tion D l mk (Q) = (lm\U(@)\lk) for the complex Wi gner D- 
functions, the following transformation rule to the real 
Wigner D-functions D^K can ^ e deduced 



E A 



Mm A/ffc D mk {Q). 



(A6) 



the Center for Theoretical Studies of ETH Zurich. 



APPENDIX A: REAL REPRESENTATION 

In the following we discuss the transformation to real 
spherical harmonics. This transformation applies not 
only to spherical harmonics, but to all quantities depend- 
ing on the angular momentum quantum numbers (Im), 
such as the Wigner D-functions and the Clebsch-Gordan 
coefficients. For the purpose of a clear notation, quan- 
tities in the complex spherical harmonics basis will be 
written with a tilde. The transformation from complex 
spherical harmonics Yi m (Q) to real spherical harmonics 
Yin(£l) is defined by 



where 



with 



Y ln {n) = ^ h nm Y lm {9), 

m= — l 



(Al) 



(A2) 



1, n>0, 
< -1, n<0, Pn 



(-1)™^, m>0, 



V 



m = 0, (A3) 



7f' m < °- 



The coefficients A nm are best represented as a matrix. 
For —2 > n, m > 2, A takes the form 



A = 



/ 


i 











i 

a/2 


\ 







i 





i 















1 















1 

a/2 





1 







V 


1 

V2 











1 

73 


) 



(A4) 



With definition ()A1|) . the 0-dependence of Yi n (6, (j>) is 
cos(n0) for n > and sin(n</>) for negative n < 0. Cre- 
ation and annihilation operators are transformed in the 



The real Wigner-D functions are explicitly given by 

% 7 ) = 2p M (3 K d L MK (9) cos(M0 + tf 7 )+ (A7) 
2a M PMP-KdM^ K {8)cos(M(f> - Kj) if u m vk = 1, 

, 7) = 2<T M (3 M pKd^ K {6) sin(M0 + tf 7 )- 
2/3 M l3-Kd^_ K (9) sin(M0 - if 7) if tr^o* = -1. 

The functions d MK (9) are the same as used for the 



Dmk( 



complex Wigner-D functions D l m 

They are tabulated in various references' 2 ^. The real 
Wigner-D functions describe the rotation of tensor oper- 



ators Tlmi such as 



4 



dk and at, which have the 



symmetries of the real spherical harmonics: 
U{Q)T LM U\<d) 



L 

E 

M'=-t 



D L M , M {Q)T LM 



(A8) 



Finally, the Clebsch-Gordan coefficients for the real 
spherical harmonics are given by 



R 



LM 



= E A M imi A 



M 2 m 2 lv Mm K&imi L 2 rn 2 > 

(A9) 

where Cf"L , m are the usual Clebsch-Gordan coeffi- 

L\m\ L 2 m 2 

cients for the complex spherical harmonics. The Clebsch- 

the propert; 

^_^L+Li+L 2 -p>LM 



Gordan coefficients Rj^m l 2 m 2 have the property 



p* LM 
rL L 1 M 1 L 2 M 2 



r>LM 

ri L 2 M 2 L 1 M 1 



L\ M± L 2 M 2 

. (A10) 

which implies that they are real if L + L\ + L2 is even 
and imaginary otherwise. Furthermore the following or- 
thogonality relation holds 



\ r>LM r/LM 
7 fir., ft/f, T.„ Mr. 



> — 5m 1 m[&m 2 m 2 - (All) 



/ j 2L L 1 M 1 L 2 M 2 L L 2 M! 1 LjMJ 
LM 



The present work is formulated in terms of the SO (3) 
symmetry where the p <S> d electron-phonon coupling in- 
volves the Clebsch-Gordan coefficients Rf 1Mo . As 
discussed in the introduction, this is equivalent to the 
t\ u ®H g electron-phonon coupling in the context of icosa- 
hedral symmetry. The icosahedral formulation involves 



14 



the icosahedral Clebsch-Gordan Coefficients tabulated in equal when using the relationships between the compo- 
Ref. l37|. The two sets of Clebsch-Gordan coefficients are nents of the IR's as given in Tab. IIVI 
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